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Abstract 



A discrete dynamical system known as Arnold's Discrete Cat Map 
(Arnold's DCM) is given by 

(x t +\,yt+i) — (x t + yt,xt + 2y t ) mod N, 

which acts on a two-dimensional square coordinate grid of size N x N . 
The defining characteristic of this map is that it has the property that 
when the N x N grid is a picture whose pixels are assigned (x, y) 
coordinates, the map scrambles the picture with each iteration. After 
a finite number of iterations, the picture is restored to its original shape 
and order. The number of iterations needed to restore the image M, 
has a mysterious dependence on N. This period, as we will find out, is 
directly related to the divisibility of the Fibonacci numbers. We will 
exploit this property to show that, for any N, an image is not dense in 
itself. In the second half of the paper, we build on the work of Chen, 
Mao, and Chui to extend the DCM to three dimensions. Finally, we 
define the generalized n dimensional DCM by introducing the idea of 
a "matrix union" . 
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1 The nature of the 2D DCM 



Arnold's Cat Mapping is a chaotic mapping on the two dimensional torus 
given by 

A(x, y) = (x + y, x + 2y) modi, 
but for our purposes, we will use the discretized version 

(x t+ i,y t +i) = (x t + yt,x t + 2y t ) mod N, 

where N is the size of the discrete grid on which the map acts. The 2D 
DCM can be written equivalently in matrix form as 

(:;;:M: ™ an - 

Example 1.1. Consider an arbitrary 3x3 image given by 

ABC 
D E F 
GUI 

The (x, ^-coordinates of the pixels in this image are 

(0,2) (1,2) (2,2)\ 
(0,1) (1,1) (2,1) mod 3. 
(0,0) (1,0) (2,0)/ 

Let's look at the orbit of an arbitrary pixel under the DCM, 

(1,1)^(2,0)^(2,2)^(1,0)^(1,1). 

The orbit of the entire image looks like 

A B C\ (B D A / D F E\ 

D E F If h a\->ia c B 

G H i) \G C E) \G I HJ 

(F A H\ (A B C\ 
-> IB I D \ \d E F . 
\G E CJ \G H I ) 

So a 3 x 3 image has period M = 4 since four iterations of the cat map were 
needed to return the image to its original state. 
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pixel dimension of image (N x N) 


iterations to restore image (period) 


300 x 300 


300 


257 x 257 


258 


183 x 183 


60 


157 x 157 


157 


150 x 150 


300 


147 x 147 


56 


124 x 124 


15 


100 x 100 


150 



Table 1: Notice, there is not an apparent correlation between the size of an 
image and its period. 



Dyson and Falk give some relationships and bounds on the period of an 
image of size N: 



M < 



= 3N for N 
= 2N for N 



< ^ for other N 



2(5 S ) ;s € N 

5 s or 6(5 S ) ; s G N 



2 An image is not dense in itself 



As in Example 1.1, consider the orbit of a single pixel as it jumps around 
an image with each iteration. An interesting question to ask is, "Does there 
exist iV such that every pixel visits every other pixel at least once?" or 
perhaps more formally, "Does there exist N such that an image is dense in 
itself?" 

To answer this question, it will be convenient to use the matrix 

A = (l l) ' Where ^ =: A<2 = (l 2) ' 

and define the Fibonacci sequence as uq = 0, u\ = 1, U2 = 1, U3 = 2,u^ = 
3, [ut+2 = ut+i + ut]. Then, the matrix for the t-th iteration of the DCM 
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is given by 




«2t+l 



) 



For a given N, the period M of the DCM is the smallest positive integer 
t, such that U2t = mod A?", and U2t-i = 1 mod N. This implies that 
u>2t+i = 'U2t+2 = 1 mod AT. Thus, the period M of the DCM for a given 
image size A" is strongly related to the divisibility of the Fibonacci numbers. 
Dyson and Falk give the following theorem and proof in [2] . The fact that 
an image is not dense in itself follows corollary. 

Theorem 2.1. For an image of size N x N, where N > 3, the period M of 
this image satisfies M < N 2 /2. 

Before we can prove this theorem, we will prove three short lemmas. 

Definition 2.2. Define fa as the least non-negative residue of itj, the iih 
Fibonacci number mod N. That is, Ui = fa mod N. 

Consider the sequence of ordered pairs (fa, fa), {fa, fa), ... , (fa, fa+i), ■ ■ ■ 
There are at most iV 2 distinct pairs. Any set of A^ 2 + 1 pairs contains some 
equal ones among them. 

Lemma 2.3. The first pair that repeats in the above sequence is (1, 1). 

Proof. Assume the opposite: that the first repeated pair is ((f) k, fa+i), where 
k > 1. So let us find a pair (fa, fa+i) where r > k in the sequence such that 
0fc = fa: ^fc+i = fa+i- From the definition of the Fibonacci numbers, 



so fa r _i = and we have that (fa-i,fa) = (fa-i, fa). But (fa-i, fa) 
occurs earlier in the sequence than (fa, fa-i)', therefore (fa, fa-i) is not the 
first pair that repeats itself. So the supposition k > 1 is wrong. Therefore, 



fa-i = fa+i ~ fa, 
fa-i = fa+i ~ fa, 



k = 1. This proves the lemma. 



□ 
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Lemma 2.4. For any positive integer N , at least one number divisible by 
N can be found among the first N 2 Fibonacci numbers. 

Proof. From the previous lemma, (1, 1} is the first pair that repeats itself. 
So (4>t, 4>t+i) = (1, 1) for some integer t such that 1 < t < N 2 + 1. Thus 



and 
But 

Therefore, 

This proves the lemma. 



(fit = 1 mod N, 

4>t+i = 1 mod N. 
ut-i = ut+i ~ ut- 
<fi t -\ = mod N. 



□ 



Lemma 2.5. For N > 2, if u t = mod N and ut+i = 1 mod N, then t 
must be even. 

Proof. The lemma is equivalent to the statement that for N > 2, if A f = 1 
mod N, then t is even. But det(A) = -1, so det(A') = (det Af = (-1)' = 1 
mod N. Hence t must be even. This proves the lemma. □ 

Theorem 2.6. The period, M of the DCM satisfies M < 

Proof. From the first and second lemmas, the second occurrence of the pat- 
tern 0, 1, 1 in the sequence <fio, <fii, (fit, (fit+i, ■■■ happens for (fij-i, (fij, <f>j+i, 
where < j — 1 < iV 2 . From the third lemma, j — 1 must be even. From 
the definition of the period, we have that 2M = j — 1. This proves the 
theorem. □ 

Corollary 2.7. An image is not dense in itself. 

Proof. Suppose there exists an iV such that w N x N image was dense in 
itself. For this to be true, it would have to be the case that the period of the 
image equals iV 2 — 1 since there are iV 2 pixels and without any iterations, the 



pixel occupies itself. But from the theorem, an arbitrary pixel in an N x N 
image has a period at most of half the total number of pixels in the image; 
a contradiction. Therefore, the number of unique places that an arbitrary 
pixel occupies in one full cycle is never more than half of the number of 
possible places. Since the choice of pixel is arbitrary, this theorem holds for 
the whole image. So an image is not dense in itself under the discrete cat 
map. This proves the corollary. □ 



3 Higher dimensional analogs 

What would a three dimensional analog to the 2D DCM look like? Certainly, 
the mapping would have to act onanA r xJVxA r cube. An iV-cube would have 
to posses a finite period, like its 2D cousin. The mapping would also have to 
posses the area-preserving and mixing dynamics of the 2D DCM. A plausible 
way to obtain an explicit form for the 3-D DCM would be to compose three 
mappings. Each mapping would fix one coordinate of x, y, z-space and have 
the lower dimensional DCM act on the remaining two coordinates. This 
way, after composition of the three "basis" maps, the entire N x N x N 
cube has been mapped in the way of the cat. 



We have 




Fixing x, we obtain the first of three 



a- 



hasis" 



maps, 




Fixing y, we obtain the second basis map, 




Fixing z, the third basis map is 



1 1 
a\ = 1 2 
\0 1 
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Composition of linear functions dictates that we multiply their respective 
matrix representations, 



A 3 = a 3 a 3 a 3 




This is the matrix for the three dimensional discrete cat map. 



3.1 The concept of "matrix union 



The systematic extrapolation of the 2D DCM to three dimensions can be 
generalized for an arbitrary positive integer dimension. We just need to 
introduce a bookkeeping device to formalize the action of fixing a coordinate 
in a basis map. 

Definition 3.1. {F^} is the set of "n x n z-frames", 

/ \ 










V 







where d% = 1 is the ith diagonal and the ith row and column except for the 
entry [di] consist entirely of O's. 



Example 3.2. 



f: 



o 
l 





N 



Ft 



Definition 3.3. The matrix union 13 is a binary function 13 : A n x F^ +1 — > 
a l n+1 which takes an n dimensional DCM matrix A n and inserts it into an 
n + 1 dimensional i-frame. The output is an n + 1 dimensional basis map 



*n+l' 



which fixes the ith coordinate of the n dimensional DCM. 
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Definition 3.4. A basis map a l n+1 is an element of the "union basis" of A n , 
= {U(A n , F^ + i)}™^ ■ This matrix fixes the ith coordinate of the n 
dimensional DCM. 

Example 3.5. 



1 0> 

Oil), 0(A 3 ,F^ 
vO 1 2j 



(1 1 1\ 

2 3 2 

10 

\3 4 4/ 



This new machinery allows for a concise iterative expression for the 
matrix for the n dimensional DCM: 

n+l 



A n+ i = I]^An,^ +1 ); n>2. 



i=l 



Example 3.6. Generate the matrix for the AD DCM. Using ([T]), we have 

(I 1 1\ 



A 













1 


1 


1 





2 


3 


2 


Vo 


3 


4 





10 
2 3 2 
\3 4 4/ 



(1 1 1\ 

2 3 2 

10 

\3 4 4/ 



4 

a 4 



(1 1 1 0\ 

2 3 2 

3 4 4 
\0 1/ 



And so 



12 3 4 
O4H4CI4O4 



/ 17 23 18 5 \ 

110 149 117 31 

257 348 274 72 

\432 585 460 122/ 



3.2 Dynamical properties of higher dimensional cat maps 



We now focus our attention on the general cat map whose t-th iteration is 
given by xt = A n *a?o, where A n is the matrix derived in the preceding sec- 
tion. In this mapping, we allow xq to be continuous. Unlike the DCM, this 
map allows for points with irrational coordinates. This way, the mapping 
has chaotic orbits. 
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Definition 3.7. (Alligood, Saur, and Yorke Def. 5.2) Let f be a map of 

W 11 , m > 1, and let {vo, Vi, V2, ...} be a bounded orbit of f. This orbit is 
chaotic if: it is not asymptotically periodic, f has no eigenvalue equal to 1, 
and at least one eigenvalue is greater than 1 in absolute value. 

For an orbit of the general cat map, an argument using folliations of the 
torus (which is beyond the scope of this paper) can be used to show that if a 
point has irrational coordinates, then its orbit is not asymptotically periodic. 
For the second two conditions, it boils down to finding the characteristic 
polynomial of A n . If (x — 1) \ XA n (A) and at least one root of XA n (A) is 
greater than 1 in absolute value, then any orbit of the general cat map whose 
coordinates are irrational is chaotic. Characteristic polynomials are given 
for the first few cases of the n dimensional general cat map below: 



Xa 2 (A) = A 2 -3A + 1, 

XA 3 (A) = -A 3 + 8A 2 -6A + 1, 

XA 4 (A) = A 4 - 562A 3 + 410A 2 - 66A + 1. 

Approximate forms for eigenvalues are given respectively as: 

0.381966, 2.61803; 

0.243019, 0.572771, 7.18421; 

0.0168808, 0.209427, 0.50397, 561.27. 

As you can see, so far there are not any eigenvalues equal to 1 and at 
least one eigenvalue is larger than 1 in absolute value in each case. These 
eigenvalues can be loosely interpreted as the amount of stretching on a per 
iteration basis of the nearby orbits under the map, where the largest value 
dominates the stretching behavior. In a sense, these numbers measure how 
chaotic a system is. 

I conjecture that as the dimension of the generalized cat map increases, 
the dominant eigenvalue of the map increases. Hence, larger dimensional 
analogs are more chaotic than their lower dimensional counterparts. Since 
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we have an iterative way using the matrix union to generate higher dimen- 
sional matrices for an arbitrary cat map which factors as a product of easily 
understood basis maps, I propose that in order to prove or disprove this, 
one would need to develop some way of studying the interaction of char- 
acteristic polynomials under the action of multiplying their corresponding 
matrices. This way, one may be able to establish a bijection between the 
natural numbers and absolute values of eigenvalues of cat maps of increasing 
dimension. 
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